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Smoothing Supernova Data to Reconstruct the Expansion 
History of the Universe and its Age 
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ABSTRACT 

We propose a non-parametric method of smoothing supernova data over redshift using 
a Gaussian kernel in order to reconstruct important cosmological quantities including 
H(z) and w{z) in a model independent manner. This method is shown to be successful 
in discriminating between different models of dark energy when the quality of data 
is commensurate with that expected from the future SuperNova Acceleration Probe 
(SNAP). We find that the Hubble parameter is especially well-determined and useful 
for this purpose. The look back time of the universe may also be determined to a 
very high degree of accuracy ( < 0.2%) in this method. By refining the method, it is 
also possible to obtain reasonable bounds on the equation of state of dark energy. We 
explore a new diagnostic of dark energy- the 'w-probe'- which can be calculated from 
the first derivative of the data. We find that this diagnostic is reconstructed extremely 
accurately for different reconstruction methods even if f2om is marginalized over. The 
u;-probe can be used to successfully distinguish between ACDM and other models of 
dark energy to a high degree of accuracy. 

Key words: cosmology: theory — cosmological parameters — statistics 



1 INTRODUCTION 

The nature of dark energy has been 
of m uch debate over th e past decade 
Sahni fc Starobinskvl 



200Ch: 



the subject 
(fo r r eview s 
Carrolll l|200ll); 



iPeebles fc Ratral (I2003D ; iPadmanabharJ (|2003h : ISahnil 



( 2004 ) 1 ). The supernova (SNe) type la data, which gave the 
first indications of the accelerated expansion of the universe, 
are expected to throw further light on this intriguing ques- 
tion as their quality steadily improves. While the number 
of SNe available to us has increased two-fold over the past 
couple of years (at present there are about 150 SNe between 
redshif ts of and 1.75 , with 10 SNe above a redshift of 
unity) (iRiess et al.|[l998l;|Perlmutter et al.|[l999l ; iKnop et all 
l2003l ; lTonrv et al.ll2003l ; iRiess et alj|2004 ). the SNe data are 
still not of a quality to firmly distinguish different models 
of dark energy. In this connection, an important role in 
our quest for a deeper understanding of the nature of dark 
energy has been played by the 'reconstruction program'. 
Commencing from t he first theoretica l exposition of the 
recons t ruction idea -IStarobinskvl (| 1998f ) ; iHuterer fc Turner! 
l|l999h : iNakamura fc Chibal (|l999l h and ISaini et all (|200Ch 



which applied it to an early supernova data set- there have 
been many attempts to reconstruct the properties of dark 
energy directly from observational data without assuming 
any particular microscopic/phenomenological model for the 
former. When using SNe data for this purpose, the main 
obstacle is the necessity to: (i) differentiate the data once 
to pass from the luminosity distance dr, to the Hubble 
parameter H(t) = a(t)/a(t) and to the effective energy 
density of dark energy eoE, (ii) differentiate the data a 
second time in order to obtain the deceleration parameter 
q = — da /a 2 , the dark energy effective pressure pde, and 
the equation of state parameter w(t) = pde/^de- Here, 
a(t) is the scale factor of a Friedmann-Robertson- Walker 
(FRW) isotropic cosmological model which we further 
assume to be spatially flat, as predicted by the simplest 
variants of the inflationary scenario of the early Universe 
and confirmed by observational CMB data. 

To get around this obstacle, some kind of smoothing 
of dr, data with respect to its argument - the redshift 
z(t) - is needed. One possible way is to parameterize the 
quantity which is of interest (H(z), w(z), etc.) by some 
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functional form containing a few free parameters and then 
determine the value of these parameters which produce the 
best fit to the data. This implies an implicit smoothing of 
cLl with a characteristic smoothing scale defined by the 
number of parameters, and with a weight depending on 
the form of parameter ization. Different parameterizations 
have been used for: dh llHuterer fc Turnerll 19991 ; ISaini et al.l 
200d; 



Alam et al 



Chib a fc Nakam ura 2000) 
20041: 



HI 



dSahni et al.l 120031; 



w(z) ( Chevallier fc Polarskil 
2002; iGerke fc Efstathioul 
2002 ; ICorasaniti fc Copelanc 
Wang fc Mukheried |2004| ; 



Alam, Salmi & Starobinskv 



2004 



Weller fc Albrecht 



a,ng & 

2004; iNesseris fc Perivelaroupolosl 



Lindoi 



Maor et al. 
20031: 



Saini. Weller fc Bridle! 
20041 ; iGond 12005 al ; 



Lazkoz. Nesseris fc Perivelaroupolosl l2005 | ') and Viz 



(Simon. 
In lHute 



Verde fc Jimenezll2005l ; iGuo, Ohta fc Zhangj|2005l ). 
Huterer fc Turner! |l999j), a polynomial expansion of the 
luminosity distance was used to reconstruct the equation 
of state. However, IWeller fc Albrecht! (|2002h showed this 
ansatz to be inadequate since it needed an arbitrarily 
large number of parameters to fit even the simplest ACDM 
equation of state. They proposed instead a polynomial 
ansatz fo r the equat i on of state which worked somewhat 
better. In lSaini et al] lj200d ) a rational Pade-type ansatz for 
d,L was proposed, which gave good results. In recent times 
there have been many more attempt s at p arame terizing 
dark energy. In IChevallier fc Polarskil (|200ll ) and iLinderl 
(2003) an ansatz of the form w = wp + w a (l — a) was 
sugges ted for the equation of state. ICorasaniti fc Copelandl 
l|2003l ) s uggested a fo u r-par ameter ansatz for the equation 
of state. ISahni et all (|2003h proposed a slightly different 
approach in which the dark energy density was expanded 
in a polyno mial ansatz. the properties of which were 
examined in dAlam et all 20041 ; lAlam. Sahni fc Starobinskv 
12004 al; lATam et al] 12004 bl) . SeelAlam et al.l (120031 ); iGone 
l|2005 bh T lBasset. Corasaniti fc Kunj l|2004h for a summary 
of different approaches to the reconstruction program and 
for a more extensive list of references. In spite of some 
ambiguity in the form of these different parameterizations, 
it is reassuring that they produce consistent results for 
the best fit curve over the range 0.1 <, z < 1 whe re we 
have su fficient amount of data (see, e.g., Fig. 10 in Gong 
l|2005 bl) ). However it is necessary to point out that the 
current SNe data are not of a quality that could allow us 
to unambiguously differentiate ACDM from evolving dark 
energy. That is why our focus in this paper will be on better 
quality data (from the SNAP experiment) which should be 
able to successfully address this important issue. 

A different, non-parametric smoothing procedure 
involves directly smoothing either d,L, or any other 
quantity defined within redshifts bins, with some 
characteristic smoothing scale. Differe nt forms of this 
appro ach have been elabor ated in IWang fc Lovelace] 



ll200ll): iHuterer fc Starkmanl J2003T ); ISainil J2003h: 



iDalv fc Diorgovskvl ll2003l. I200I) 



2005); Espana-Bonet fc Ruiz-Lapuentd 



Wang fc Tegmarkl 
20051 ). One of the 



advantages of this approach is that the dependence of the 
results on the size of the smoothing scale becomes explicit. 
We emphasize again that the present consensus seems to 
be that, while the cosmological constant remains a good 
fit to the data, more exotic models of dark energy are 
by no means ruled out (though their diversity has been 



significantly narrowed already). Thus, until the quality 
of data improves dramatically, the final judgment on the 
nature of dark energy cannot yet be pronounced. 

In this paper, we develop a new reconstruction method 
which formally belongs to the second category, and which 
is complementary to the approach of fitting a parametric 
ansatz to the dark energy density or the equation of state. 
Most of the papers using the non-parametric approach cited 
above exploited a kind of top-hat smoothing in redshift 
space. Instead, we follow a procedure which is well known 
and frequently used in the analysis of large-sc ale structure 
l|Coles fc Lucchinlll995l : iMartinez fc Saaril2002ft : namely, we 
attempt to smooth noisy data directly using a Gaussian 
smoothing function. Then, from the smoothed data, we cal- 
culate different cosmological functions and, thus, extract 
information about dark energy. This method allows us to 
avoid additional noise due to sharp borders between bins. 
Furthermore, since our method does not assume any defi- 
nite parametric representation of dark energy, it does not 
bias results towards any particular model. We therefore ex- 
pect this method to give us model-independent estimates of 
cosmological functions, in particular, the Hubble parameter 
H(z) = a(t)/a(t). On the basis of data expected from the 
SNAP satellite mission, we show that the Gaussian smooth- 
ing ansatz proposed in this paper can successfully distin- 
guish between rival cosmological models and help shed light 
on the nature of dark energy. 



2 METHODOLOGY 

It is useful to recall that, in the context of structure for- 
mation, it is often advantageous to obtain a smoothed den- 
sity field S s (x) from a fluctuating 'raw' density field, 5(x'), 
using a low pass filter F having a characteristic scale Rf 
|Coles fc Lucchinlll995l ) 



(x,7?/)= j 8(-K)F(\x.-x.'\;Rf) dx' 



(1) 



Commonly used filters include: (i) the 'top-hat' filter, which 
has a sharp cutoff Fth oc 0(1 — |x — x'|/_Rth), where 
is the Heaviside step function (0(z) = for z ^ 0, 
0(2) = 1 for z > 0) and (ii) the Gaussian filter Fq oc 
exp(— |x — x'| 2 /2Rq). For our purpose, we shall find it use- 
ful to apply a variant of the Gaussian filter to reconstruct 
the properties of dark energy from supernova data. In other 
words, we apply Gaussian smoothing to supernova data 
(which is of the form {lndi(zi), Zi}) in order to extract in- 
formation about important cosmological parameters such as 
H(z) and w(z). The smoothing algorithm calculates the lu- 
minosity distance at any arbitrary redshift z to be 

\nd L (z,A) s = ln d L (z) 9 + N(z)Y^ ?nd L (z t ) -In d L (z t ) 9 ] 



x exp 



i* 2 m) 



2A 2 



N(z)-i=J2 



exp 



i» 2 mi 

2A 2 



Here, lnd_L(z,A) s is the smoothed luminosity distance at 
any redshift z which depends on luminosity distances of each 
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Table 1. Expected number of supernovae per redshift bin from the SNAP experiment 



Az 0.1-0.2 0.2-0.3 0.3-0.4 0.4-0.5 0.5-0.6 0.6-0.7 0.7-0.8 0.8-0.9 
N 35 64 95 124 150 171 183 179 



Az 0.9-1.0 1.0-1.1 1.1-1.2 1.2-1.3 1.3-1.4 1.4-1.5 1.5-1.6 1.6-1.7 
N 170 155 142 130 119 107 94 80 



SNe event with the redshift Zj, and N(z) is a normaliza- 
tion parameter. Note that the form of the kernel bears re- 
semblance to the lognormal distribution (such distributions 
find applic ation in the study of cosmological density per- 
turbations, ISahni fc Colesl (|l995h V The quantity lnd L (z) 9 
represents a guessed background model which we subtract 
from the data before smoothing it. This approach allows 
us to smooth noise only, and not the luminosity distance. 
After noise smoothing, we add back the guess model to re- 
cover the luminosity distance. This procedure is helpful in 
reducing noise in the results. Since we do not know which 
background model to subtract, we may take a reasonable 
guess that the data should be close to ACDM and use 
cIl{z) 9 = di(z) ACDM as a first approximation and then use 
a boot-strapping method to find successively better guess 
models. We shall discuss this issue in greater detail in the 
section [3] Having obtained the smoothed luminosity dis- 
tance, we differentiate once to obtain the Hubble param- 
eter H(z) and twice to obtain the equation of state of dark 
energy w(z), using the formula 



H(z) = 



d 

dz 



d L (z) 
l + z 



[2(l + z)/3] H'/H-l 

i - (# /-ff) 2 n 0m (l + «) 3 



(3) 



(4) 



The results will clearly depend upon the value of the scale 
A in ([2]). A large value of A produces a smooth result, 
but the accuracy of reconstruction worsens, while a small 
A gives a more accurate, but noisy result. Note that, for 
\z — Zi\ <C 1, the exponent in Eq. ([2| reduces to the form 
— (z — Zi) 2 /2A 2 (l+z) 2 . Thus, the effective Gaussian smooth- 
ing scale for this algorithm is A(l + z). We expect to obtain 
an optimum value of A for which both smoothness and ac- 
curacy are reasonable. 

The Hubble parameter can also be used to obtained the 
weighted average of w 



1+w 



51n(l + z) 



:i+w(z)) 



dz 



1 S InpDE 



l + z 3 8 ln(l + z) 



.(5) 



Pde is the dark energy density pDE = Pde/poc (where 
poc = 3Hq /8-kG). We shall show in the section [S] that w, 
which we call the w-probe, acts as an excellent diagnostic of 
dark energy, and can differentiate between different models 
of dark energy with greater accuracy than the equation of 
state. 

To check our method, we use data simulated accord- 
ing to the SuperNova Acceleration Probe (SNAP) experi- 
ment. This space-based mission is expected to observe close 
to 6000 supernovae, of which about 2000 supernovae can be 



used for cosmological purposes (|Aldering et al.ll2004l ). We 
propose to use a distribution of 199 8 supernovae betwee n 
redshifts of 0.1 and 1.7 obtained from lAldering et all l|2004h . 
This distribution of 1998 supernovae is shown in Table Q] 
Although SNAP will not be measuring supernovae at red- 
shifts below z — 0.1, it is not unreasonable to assume that, 
by the time SNAP comes up, we can expect high quality 
data at low redshifts from other supernova surveys such as 
the Nearby SN Factory 0. Hence, in the low redshift region 
z < 0.1, we add 25 more supernovae of equivalent errors to 
the SNAP distribution, so that our data sample now con- 
sists of 2023 supernovae . Using this distribution of data, 
we check whether the method is successful in reconstruct- 
ing different cosmological parameters, and also if it can help 
discriminate different models of dark energy. 

We simulate 1000 realizations of data using the SNAP 
distribution with the error in the luminosity distance 
given by a\ n d L = 0.07 - the expected error for SNAP. 
We also consider the possible effect of weak-lensing on 
high redshift supernovae by adding an uncertainty of 
aiens(z) fts 0.46( 0031 1 + 0.08687z - 0.00950z 2 ) (as in 
IWang fc Tegmarkl |2005.)). Initially, we use a simple model 
of dark energy when simulating data - an evolving model of 
dark energy with w = —a/ao = —1/(1 + z) and fiom = 0.3. 
It will clearly be of interest to see whether this model can 
be reconstructed accurately and discriminated from ACDM 
using this method. From the SNAP distribution, we obtain 
smoothed data at 2000 points taken uniformly between the 
minimum and maximum of the distributions used. Once we 
are assured of the efficacy of our method, we shall also at- 
tempt to reconstruct other models of dark energy. Among 
these, one is the standard cosmological constant (ACDM) 
model with w = — 1. The other is a model with a constant 
equation of state, w = —0.5. Such models with constant 
equatio n of state are kn own as quiessence models of dark 
energy |Alam et alj|2003h and we shall refer to this model as 
the "quiessence model" throughout the paper. These three 
models are complementary to each other. For the ACDM 
model, the equation of state is constant at w = — 1, w re- 
mains constant at —0.5 for the quiessence model and for the 
evolving model, w(z) varies rapidly, increasing in value from 
Wo = — 1 at the present epoch to w — at high redshifts. 



http: / /snfactory.lbl.gov 
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A = 0.24 
Fiducial Model: w = -1/(1 + z) 




Figure 1. The smoothing scheme of equation J2j is used to determine H(z) and w(z) from 1000 realizations of the SNAP dataset. The 
smoothing scale is A = 0.24. The dashed line in each panel represents the fiducial w = —1/(1 + z) 'metamorphosis' model while the solid 
lines show the mean Hubble parameter (left), the mean equation of state (right), and \a limits around these quantities. The dotted line 
in both panels is ACDM. Note that the mean Hubble parameter is reconstructed so accurately that the fiducial model (dashed line) is 
not visible in the left panel. 



3 RESULTS 

In this section we show the results obtained when our 
smoothing scheme is applied to data expected from the 
SNAP experiment. The first issue we need to consider is that 
of the guess model. As mentioned earlier, the guess model in 
equation ([2]) is arbitrary. Using a guess model will naturally 
cause the results to be somewhat biased towards the guess 
model at low and high redshifts where there is paucity of 
data. Therefore we use an iterative method to estimate the 
guess model from an initial guess. 
Iterative process to obtain Guess model 

To estimate the guess model for our smoothing scheme, 
we use the following iterative method. We start with a sim- 
ple cosmological model, such as ACDM, as our initial guess 
model- \ndf = lnd£ CDM . The result obtained from this 
analysis, lnd^, is expected to be closer to the real model 
than the initial guess. We now use this result as our next 
guess model- In d 9 L = In d\ and obtain the next result In d 2 L . 
With each iteration, we expect the guess model to become 
more accurate, thus giving a result that is less and less biased 
towards the initial guess model used. A few points about the 
iterative method should be noted here. 

• Using different models for the initial guess does not af- 
fect the final result provided the process is iterated several 
times. For example, if we use a w = —1/(1 + z) 'metamor- 
phosis' model to simulate the data and use either ACDM 
or the w = —0.5 quiessence model as our initial guess, the 
results for the two cases converge by > 5 iterations. 

• Using a very small value of A will result in a accurate 
but noisy guess model, therefore after a few iterations, the 
result will become too noisy to be of any use. Therefore, 



we should use a large A for this process in order to obtain 
smoother results. 

• The bias of the final result will decrease with each it- 
eration, since with each iteration we get closer to the true 
model. The bias decreases non-linearly with the number of 
iterations M. Generally, after about 10 iterations, for mod- 
erate values of A, the bias is acceptably small. Beyond this, 
the bias still decreases with the number of iterations but the 
decrease is negligible while the process takes more time and 
results in larger errors on the parameters. 

• It is important to choose a value of A which gives a 
small value of bias and also reasonably small errors on the 
derived cosmological parameters. To estimate the value of 
A in we consider the following relation between the 
reconstructed results, quality and quantity of the data and 
the smoothing parameters. One can show t hat the relative 
error bars on H(z) scale as jT^gmark|[2002) 



H A/1/2A3/2 ' v ' 

where N is the total number of supernovae (for approxi- 
mately uniform distribution of supernovae over the redshift 
range) and a is the noise of the data. From the above equa- 
tion we see that a larger number of supernovae or larger 
width of smoothing, A, will decrease the error bars on re- 
constructed H, but as we shall show in appendix [A] the 
bias of the method is approximately related to A 2 . This im- 
plies that, by increasing A we will also increase the bias of 
the results. We attempt to estimate A such that the error 
bars on H be of the same order as a, which is a reasonable 
expectation. 

If we consider a single iteration of our method, then for 
N ~ 2000 we get A ~ N~ 1/3 ~ 0.08. However, with each 
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A = 0.24 

Fiducial Model: w = -1/(1+ z) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



z 

Figure 2. The smoothing scheme of equation (f2j) is used to de- 
termine the look-back time of the universe, T(z) = t(0) — t(z), 
from 1000 realizations of the SNAP dataset for a w = —1/(1 + z) 
'metamorphosis' model. The smoothing scale is A = 0.24. The 
solid lines show the mean look-back time and the la limits around 
it. The look-back time for the fiducial model matches exactly with 
the mean for the smoothing scheme. The dotted line shows the 
ACDM model. 

iteration, the errors on the parameters will increase. There- 
fore using this value of A when we use an iterative process 
to find the guess model will result in such large errors on 
the cosmological parameters as to render the reconstruction 
exercise meaningless. It shall be shown in Appendix lAl that 
at the M-th iteration, the error on lndi, will be approxi- 
mately 8 At (in <1l) — VmSoQucIl)- The error on lndz, scales 
as 1/A. We would like the errors after M iterations to be 
commensurate with the optimum errors obtained for a single 
iteration, Ao, so we require A op ti ma i — vMAo. Therefore, if 
we wish to stop the boot-strapping after 10 iterations, then 
Aoptimai — 3Ao — 0.24. This is the optimal value of A we 
shall use for best results for our smoothing procedure. 

Considering all these factors, we use a smoothing scale 
A = 0.24 for the smoothing procedure of Eq ((2)1 with a itera- 
tive method for finding the guess model (with ACDM as the 
initial guess). The boot-strapping is stopped after 10 itera- 
tions. We will see that the results reconstructed using these 
parameters do not contain noticeable bias and the errors on 
the parameters are also satisfactory. 

Figure [1] shows the reconstructed H(z) and w(z) with 
la errors for the w = —1/(1 + z) evolving model of dark 
energy. From this figure we can see that the Hubble param- 
eter is reconstructed quite accurately and can successfully 
be used to differentiate the model from ACDM. The equa- 
tion of state, however, is somewhat noisier. There is also a 
slight bias in the equation of state at low and high redshifts. 
Since the w — —1/(1 + z) model has an equation of state 
which is very close to w = — 1 at low redshifts, we see that 
w(z) cannot discriminate ACDM from the fiducial model at 
z < 0.2 at the la confidence level. 



Age of the Universe 

We may also use this smoothing scheme to calculate 
other cosmological parameters of interest such as the age of 
the universe at a redshift z: 
r°° j j 

In this case, since data is available only upto redshifts of 
z ~ 1.7, it will not be possible to calculate the age of the 
universe. Instead, we calculate the look-back time at each 
redshift- 

T(*)=t(z = 0)-t(z) = H ^J* {l+ d ; mz/y (8) 

Figure[5]shows the reconstructed T(z) with la errors for 
the w = —1/(1 + 2;) 'metamorphosis' model using the SNAP 
distribution. For this model the current age of the universe 
is about 13 Gyrs and the look-back time at z ~ 1.7 is about 
9 Gyrs for a Hubble parameter of -Ho = 70 km/s/Mpc. We 
see that the look-back time is reconstructed extremely ac- 
curately. Using this method we may predict this parameter 
with a high degree of success and distinguish between the 
fiducial look-back time and that for ACDM even at the 10cr 
confidence level. Indeed any cosmological parameter which 
can be obtained by integrating the Hubble parameter will be 
reconstructed without problem, since integrating involves a 
further smoothing of the results. 

Looking at these results, we draw the conclusion that 
the method of smoothing supernova data can be expected to 
work quite well for future SNAP data as far as the Hubble 
parameter is concerned. Using this method, we may recon- 
struct the Hubble parameter and therefore the expansion 
history of the universe accurately. We find that the method 
is very efficient in reproducing H(z) to an accuracy of <, 2% 
within the redshift interval < z < 1, and to <, 4% at 
z ~ 1.7, as demonstrated in figure [T] Furthermore, using 
the Hubble parameter, one may expect to discriminate be- 
tween different families of models such as the metamorphosis 
model w — — l/(l + z) and ACDM. This method also repro- 
duces very accurately the look-back time for a given model, 
as seen in fig [5] It reconstructs the look-back time to an 
accuracy of < 0.2% at z ~ 1.7. 



4 REDUCING NOISE THROUGH DOUBLE 
SMOOTHING 

As we saw in the preceding section, the method of smoothing 
supernova data to extract information on cosmological pa- 
rameters works very well if we employ the first derivative of 
the data to reconstruct the Hubble parameter. It also works 
reasonably for the second derivative, which is used to deter- 
mine w(z), but the errors on w(z) are somewhat large. In 
this section, we examine a possible way in which the equa- 
tion of state may be extracted from the data to give slightly 
better results. 

The noise in each parameter translates into larger noise 
levels on its successive derivatives. We have seen earlier that, 
using the smoothing scheme ([2]), one can obtain H(z) from 
the smoothed cLl(z) fairly successfully. However, small noises 
in H(z) propagate into larger noises in w(z). Therefore, it 
is logical to assume that if H(z) were smoother, the resul- 
tant w(z) might also have smaller errors. So, we attempt to 
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A = 0.24, Double Smoothing 
Fiducial Model: w = -1/(1 + z) 





Z Z 

Figure 3. The double smoothing scheme of equations J2) and (0 has been used to obtain H (z) and w[z) from 1000 realizations of the 
SNAP dataset. The smoothing scale is A = 0.24. The dashed line in each panel represents the fiducial w = — 1/(1 + z) 'metamorphosis' 
model while the solid lines represent the mean and la limits around it. The dotted line in both panels is ACDM. In the left panel H(z) 
for the fiducial model matches exactly with the mean for the smoothing scheme. 



A = 0.24, Double Smoothing 
Fiducial Model: w = — 1 




d 




Z Z 

Figure 4. The double smoothing scheme of equations (|2j and (O has been used to obtain H(z) and w(z) from 1000 realizations of the 
SNAP dataset. The smoothing scale is A = 0.24. The dashed line in each panel represents the fiducial ACDM model with w = —I while 
the solid lines represent the mean and la limits around it. In the left panel H(z) for the fiducial model matches exactly with the mean 
for the smoothing scheme. 



smooth H(z) a second time after obtaining it from cLl(z). 
The procedure in this method is as follows - first, we smooth 
noisy data lndi(z) to obtain lndi(2) s using equation ([2]). 
We differentiate this to find H(z) s using equation We 
then further smooth this Hubble parameter by using the 



same smoothing scheme at the new redshifts 
H(z, A) 52 = H{zf + N(z) J2 [Hi^T - H&Y] 

i» 2 m) 



x exp 



2A 2 



(9) 
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A = 0.24, Double Smoothing 
Fiducial Model: w = -0.5 



0.5 



1 .5 




Z Z 

Figure 5. The double smoothing scheme of equations (|2j and (0 has been used to obtain H(z) and w(z) from 1000 realizations of the 
SNAP dataset. The smoothing scale is A = 0.24. The dashed line in each panel represents the fiducial quiessence model with w = —0.5 
while the solid lines represent the mean and la limits around it. The dotted line is ACDM. 



N(z)- 1 = ^exp 



In 2 



' l+Zi ' 
, 1 + 2 J 



2A 2 



We then use this H(z, A) S2 to obtain w(z) using equa- 
tion Q. This has the advantage of making w(z) less noisy 
than before, while using the same number of parameters. 
However, repeated smoothing can also result in the loss of 
information. 

The result for the SNAP distribution using this double 
smoothing scheme for the w — —1/(1 + z) model is shown 
in figure [3] We use A = 0.24 for smoothing both lndzT-z) 
and H(z). Comparing with figure [T] we find that there is an 
improvement in the reconstruction of H(z) as well as w(z). 
Thus, errors on the Hubble parameter decrease slightly and 
errors on w(z) also become somewhat smaller. 

We now explore this scheme further for other models of 
dark energy. We first consider a w = — 1 ACDM model. In 
figure I?] we show the results for this model. We find that 
the Hubble parameter accurately reconstructed and even w 
is well reconstructed, with a little bias at high redshift. The 
next model we reconstruct is a w — —0.5 quiessence model. 
The results for double smoothing are shown in fig [5] There 
is a little bias for this model at the low redshifts, although 
it is still well within the error bars. 

We note that in all three cases, a slight bias is notice- 
able at low or high redshifts. This is primarily due to edge 
effects- since at low (high) redshift, any particular point will 
have less (more) number of supernovae to the left than to 
the right. Even by estimating the guess model through an 
iterative process, it is difficult to completely get rid of this 
effect. In order to get rid of this effect, we would require to 
use much larger number of iterations for the guess model, 
but this would result in very large errors on the parameters. 
However, this bias is so small as to be negligible and cannot 
affect the results in any way. 



Looking at these three figures, we can draw the follow- 
ing conclusions. The Hubble parameter is quite well recon- 
structed by the method of double smoothing in all three 
cases while the errors on the equation of state also decrease. 
At low and high redshifts, a very slight bias persists. Despite 
this, the equation of state is reconstructed quite accurately. 
Also, since the average error in w(z) is somewhat less than 
that in the single smoothing scheme (figure [TJ, the equation 
of state may be used with better success in discriminating 
different models of dark energy using the double smoothing 
procedure. 



5 THE ro-PROBE 

In this section we explore the possibility of extracting infor- 
mation about the equation of state from the reconstructed 
Hubble parameter by considering a weighted average of the 
equation of state, which we call the w-probe. An important 
advantage of this approach is that there is no need to go 
to the second derivative of the luminosity distance for in- 
formation on the equation of state. Instea d, we consider the 
weighted average of the equation of state (|Alam et alj|2004h 



1 + w 



1 



<51n(l + z) 



(1 + «;(*)) 



dz 

m 



(10) 



which can be directly expressed in terms of the differ- 
ence in dark energy density pde = Pde/poc (where po? — 
over a range of redshift as 



i + w(zi, z 2 y 



1 S In pde 
3 5 ln(l + z) 



■4* 



h 2 (zi) -n 0m (n-«i) 3 



H 2 (z 2 )-Slom(l + Z 2 ) 3 



'In 



l + Zl 
1+22 



(11) 
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Table 2. The reconstructed ui-probe (w), determined using Eq. JTTJ (and its lc error) is listed for 1000 realizations of SNAP data. 
Three fiducial models are used : w = —1/(1 + z), w = —1 (ACDM) and to = —0.5. We deploy the method of double smoothing with 
A = 0.24 to determine w. 





w = -1/(1 - 




w = —1 




w = -0.5 




Az 


w 


^cxact 


w 


Inexact 


w 


^cxact 


- 0.414 


-0.839 ±0.019 


-0.845 


-1.001 ±0.017 


-1.0 


-0.489 ± 0.025 


-0.5 


0.414 - 1 


-0.595 ±0.033 


-0.598 


-1.009 ± 0.038 


-1.0 


-0.506 ± 0.039 


-0.5 


1 - 1.7 


-0.471 ±0.069 


-0.432 


-1.017 ±0.087 


-1.0 


-0.493 ± 0.075 


-0.5 



where 5 denotes the total change of a variable between in- 
tegration limits. Thus, even if the equation of state is noisy, 
the w parameter may be obtained accurately provided the 
Hubble parameter is well constructed. 

The parameter w has the interesting property that for 
the concordance ACDM model, it equals —1 in all redshift 
ranges while for other models of dark energy it is non-zero. 
For (non-ACDM) models with constant equation of state, 
this parameter is a constant (but not equal to —1), while 
for models with variable equation of state, it varies with 
redshift. The fact that ACDM is a fixed point for this quan- 
tity may be utilized to differentiate between the concordance 
ACDM model and other models of dark energy. Therefore 
the parameter w may be used as a new diagnostic of dark en- 
ergy which acts as a discriminator between ACDM and other 
models of dark energy. We call this diagnostic the w-probe. 

We now calculate the w-probe for the three models 
described above using the method of double smoothing. In 
tableEl we show the values of w obtained in different redshift 
ranges after applying double smoothing on SNAP-like data. 
The ranges of integration are taken to be approximately 
equally spaced in ln(l + z). Two points of interest should 
be noted here: a) w is very close to TO CX act at all redshifts 
for all three models of dark energy; and b) as expected, 
this parameter is good at distinguishing between ACDM and 
other dark energy models. 

In the above analysis, we have assumed that the mat- 
ter density is known exactly, f2om = 0.3. Studies of large 
scale structure and CMB have resulted in very tight bounds 
on the matter density, but still som e uncertainty rema ins 
regarding its true value. As noted in iMaor et al.l (|2002l ). a 
small uncertainty in the value of f2om may affect the recon- 
struction exercise quite dramatically. The Hubble parameter 
is not affected to a very high degree by the value of mat- 
ter density, because it can be calculated directly as the first 
derivative of the luminosity distance, which is the measured 
quantity. However, when calculating the equation of state of 
dark energy, the value of f2om appears in the denominator of 
the expression @, hence any uncertainty in ilom is bound 
to affect the reconstructed w(z). This is illustrated in figure 
[6] where the equation of state is reconstructed for a fiducial 
ACDM model with Qom = 0.3 by assuming that the matter 
density has been incorrectly determined to be f2om = 0.2, 
and this incorrect value is then used to determine w(z). For 
this purpose we consider a complementary reconstruction 
ansatz that uses a polynomial fit to the dark energy density 
jSahni et alj|2003h 




in 

i i i i i i i i i i i i i i i i 

0.5 1 1.5 

Z 

Figure 6. The reconstructed equation of state w(z) is shown for 
1000 realizations of a Q m = 0.3, w = -1 ACDM model. We 
assume an incorrect value for the matter density, f2o m = 0.2 in 
the reconstruction exercise performed using (J4j and the ansatz 
112> . This is done to study the effect of the observed uncertainty 
in Qom on the equation of state. The dashed line represents the 
fiducial ACDM model with w = —1.0. The solid lines represent 
the mean w(z) and the la limits around it. 



H 2 {z) =Hi[Q 0m {l + z) :i +A + A 1 (l + z)+A2{l + z) 2 ] ,(12) 

where Aq — 1 — fiom — A\ — A2 for a fiat universe. This ansatz 
is known to give accurate results for the dark energy d ensity 
( Ala m et al, .1 120041 :1X1 am. Sahni fc Starobinskvll2004 ah . Fig- 
ure [6] clearly shows that an incorrect value of Qom gives rise 
to an erroneously evolving equation of state of dark energy 
whereas, in fact, the correct EOS remains fixed at w = — 1 
and does not evolve. 

One of the main results of this paper is that, although 
the equation of state w(z) may be reconstructed badly if 
f2om is not known accurately, the uncertainty in f2o m does 
not have such a strong effect on the reconstruction of the 
w-probe w. This is because w in equation is a difference 
of two terms, both involving flom- As a result, uncertainty in 
fiom does not affect w as much as it affects w(z). Therefore, 
even when Oo m is not known to a high degree of accuracy, 
the ui-probe may still be reconstructed fairly accurately. 
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z z 

Figure 7. The ui-probe is reconstructed for the unevolving ACDM model with w = — 1 (left panel) and an evolving DE model with 
w = —1/(1 + z) (right panel). 1000 realizations of SNAP-like data have been used. The thick dashed line in both panels indicates the 
exact value of to for the fiducial model, the dark grey boxes in each panel indicate the la confidence levels on £D reconstructed for the 
two models using the double smoothing scheme with A = 0.24 and marginalising over Qom = 0.3 ± 0.07. This figure illustrates that the 
ui-probe works remarkably well for both ACDM (left panel) and for evolving DE (right panel). The details for this figure are given in 
Table [3] 



We now demonstrate this by showing the results ob- 
tained using our smoothing scheme after marginalising over 
the matter density. We simulate SNAP like data for two 
models: (a) ACDM and (b) aic = —1/(1 + 2) 'metamorpho- 
sis' model. When applying the smoothing scheme, we as- 
sume that S7om follows a Gaussian probability distribution 
with mean ttom = 0.3 and variance a = 0.07 (the error be- 
ing commensurate to that expect ed from the current CMB 
and Large Scale Structure data (|Percival et al.ll200ll )). In 
figure [7] and table O we show the results for the «;-probe 
calculated for the two models. We find that the ui-probe 
(w) is determined to a high degree of accuracy for both the 
models even when we marginalize over f2om ! The value of 
w for the ACDM model is approximately equal to —1, while 
that for the metamorphosis model shows clear signature of 
evolution. Thus, even if the matter density of the universe 
is known uncertainly, this uncertainty does not affect the 
accuracy of the reconstructed w-probe significantly. This is 
a powerful result since it indicates that unlike the equation 
of state, the «;-probe is not overtly sensitive to the value of 
flom for SNAP-quality data. 

From the above results, we see that the «;-probe is very 
effective as a diagnostic of dark energy, especially in differ- 
entiating between ACDM and other models of dark energy. 
We summarise some important properties of the w-probe 
below: 

(i) w(zi, Z2) is determined from the first derivative of the 
luminosity distance. Its reconstructed value is therefore less 
noisy than the equation of state w(z) (which is determined 
after differentiating cIl(z) twice; compare Q and (| 1 1 1) ). 

(ii) w(z\, Z2) = —1 uniquely for concordance cosmology 
(ACDM). For all other dark energy models w ^ — 1. This 
remains true when w is marginalized over fiom- 

(iii) w is robust to small uncertainties in the value of 



Table 3. The reconstructed ui-probe id (Eq. mil l over specified 
redshift ranges (and its la error) is shown for 1000 realizations of 
SNAP data. Two fiducial models are used : the w = —1/(1 + z) 
'metamorphosis' model and w = —1 (ACDM). We deploy the 
method of double smoothing with A = 0.24 and marginalize over 
Q m = 0.3 ± 0.07. 





w = -1/(1 - 




w = —1 




Az 


w 


Inexact 


w 


^cxact 


- 0.414 


-0.837 ±0.025 


-0.845 


-1.003 ±0.021 


-1.0 


0.414 - 1 


-0.618 ±0.042 


-0.598 


-1.018 ±0.052 


-1.0 


1 - 1.7 


-0.461 ±0.127 


-0.432 


-1.051 ±0.147 


-1.0 



the matter density. As we saw earlier this uncertainty can 
induce large errors in deter minations of the c osmic equa- 
tion of state w(z), see also iMaor et al.l (|2002l ). The weak 
dependence of w on the value of fiom m t ne range cur- 
rently favored by observations 0.2 ^ flo m ^ 0.4 implies 
that the ui-probe can cope very effectively with the ex- 
isting uncertainty in the value of the matter density for 
SNAP-quality data. Furthermore, since w is constructed di- 
rectly from Pde, any method which determines either the 
dark energy density or the Hubble parameter from obser- 
vations can be used to also determine w. Note that sev- 
eral excellent methods for deter mining pde and H(z) have 
been suggested in the l iterature |b alv & Diorgovsky 120031 
| 2004l:lAlam et al.ll20ol lAlam. Sahni fc Starobinskvll2004 al : 
I Wang fc Mukherieell2004l : IWang fc Tegmarkll2005l ), and any 
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of these could be used to great advantage in determining the 
w-probe. 

Thus we expect that the w-probe may be used as a 
handy diagnostic for dark energy, especially in discriminat- 
ing between ACDM and other models of dark energy, for 
SNAP like datasets. Its efficacy lies in the fact that it is 
not very sensitive to both the value of the present matter 
density and also the reconstruction method used. 



6 COSMOLOGICAL RECONSTRUCTION 

APPLIED TO OTHER PHYSICAL MODELS 
OF DARK ENERGY 

In this section we draw the readers attention to the dangers 
encountered during cosmological reconstruction of atypical 
dark energy models. There are currently two plausible 
ways of making the expansion of the universe accelerate 
at late times. The first approach depends on changing 
the matter sector of the Einstein equations. Examples of 
this approach are the quintessence fields. A completely 
different approach has shown that it is possible to obtain 
an accelerating universe th rough modifying the gravity 
sector (se e, for instance , iDeffayet. Dyali fc Gabadadze 
iFreese fc Lewisl (I2002TI; ISahni fc Shtanov 
!): ICarroll. Hoffman fc Trodden! (l2003h ; 



ICapozziello. Carloni fc Troisil 



1 2003ft : iDolgov fc Kawasaki! 



.Noiiri fc Odintsovl 
ISahnil (2005ft and 



references therein). In these models, dark energy should not 
be treated as a fluid or a field. Instead, it may be better 
dubbed as 'geometric dark energy'. Indeed the DGP model 
can cause the universe to accelerate even in the absence 
of a physic a l dar k energy comp onent. As pointed out in 
I Alam et all (120031 ); ISahnil i|2005ft . the equation of state is 
not a fundamental quantity for geometric dark energy. E.g., 
using w(z) in the reconstruction of such models may result 
in very strange results, including, for instance, singularities 
in the equation of state. 

As an exampl e, we consider the branew orld dark energy 
model proposed in lSahni fc Shtanov! (|2003T ) described by the 
following set of equations for a flat universe : 



Q Qm (l + z) + fi CT + 



(13) 



2fi; - 2Vfh v /fi 0m (l + ^) 3 + Ox + fi; + f2; b 

= 1 - ci 0m + 2 v / n l (i + n h ) , (14) 



A very simple model which has a well-behaved q(z) but sin- 
gular ui(z) is a model which has, in addition to the cosmological 
constant, a second dark energy component disguised as a spatial 
curvature term- H 2 (z)/H$ = S7 m(l + z) 3 + + -z) 2 + ^A- 

If we assume that Com = 0.3, Qx = —0.05, £1a = 0.75, then w(z) 
becomes singular when Qx (1 + z) 2 = 0, i.e., at z ~ 2.8. Al- 

though this property of w(z) can be easily understood physically 
and rests in the fact that it is an 'effective' equation of state for 
the combination of DE fluids, nevertheless any reasonable param- 
eterization of w(z) will clearly experience difficulty in reproducing 
this behavior. An effective equation of state with a similar 'pole- 
like' divergence is frequently encountered in braneworld models 
of dark energy (Sahni & Shtanov 2003, 2005J) as well as in holo- 
graphic models jLinderll2004l) . 



A = 0.24, Double Smoothing 
Fiducial Model: Braneworld : fiom = 0-3, 0; = 1, Sli b = 

LT> 



LO 

o 




z 

Figure 8. The reconstructed Hubble parameter H(z) for the 
braneworld model H13I I for 1000 realizations of the SNAP dataset. 
We use the double smoothing scheme of equations (|2j and 
with A = 0.24. The dashed line represents the fiducial f!om = 
0.3, fi; = 1, f2; b = braneworld model while the solid lines repre- 
sent the mean and la limits around it. We note that the Hubble 
parameter is very well reconstructed for the braneworld model. 



where the densities Q are defined as 



Pot> 



3m 2 H 2 ' 



3m? Hi ' 



Qi = wi^=-^k ,(15) 



l c = m 2 /M 3 being a new length scale (m and M refer re- 
spectively to the four and five dimensional Planck masses), 
lb the bulk cosmological constant and a the brane tension. 
In this section we have used % = c = 1. On short length 
scales r <C l c and at early times, one recovers general rel- 
ativity, whereas on large length scales r > i t and at late 
times brane-related effects become important and may lead 
to the acceleration of the universe. The 'effective' equation 
of state for this braneworld model is given by 

P= S (1 "" m(z)) ' P= ^G {q{z) ~ 1,2) (16) 
p q(z) - 1/2 

Weff = - = -Tf. „ , ■ (17) 

p 3(1-S2 m (z)) 

It is obvious that the effective equation of state in this 
braneworld model may become singular if Q m (z) = f2om(l + 
z ) 3 Hq I H 2 ( z ) becomes unity. This does not signal any inher- 
ent pathologies in the model however. We should remember 
that the acceleration of the universe in this model is due to 
modification of the expansion of the universe at late times 
due to extra-dimensional effects. Hence it is not very appro- 
priate to describe dark energy by an equation of state for 
such a model. However, it would be interesting to see if the 
singularity in the effective w for this model can be recovered 
by our smoothing method. 

We attempt to reconstruct an Slom = 0.3, Q,i = 1, Qi b = 
branewo rld model which is a g ood fit to the current super- 
nova data lAlam fc Sahnil |2002ft . We simulate data accord- 
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Figure 9. Reconstructed equation of state for the braneworld model <1 1 3 tt for 1000 realizations of the SNAP dataset. The left panel shows 
results for the double smoothing scheme of equations |(2j and (0 with A = 0.24 while the right panel uses the ansatz il8l to reconstruct 
the equation of state. The dashed line in each panel represents the fiducial f2om = 0.3, Qi = 1, Qi b = braneworld model while the solid 
lines represent the mean and lcr limits around it. 



ing to SNAP and obtain results for the double smoothing 
method with A = 0.24. In figure |8j we show the recon- 
structed Hubble parameter for this reconstruction. We see 
that the Hubble parameter is very well reconstructed and 
shows no pathological behavior. 

We now obtain the equation of state of dark energy 
for this model. For this purpose, we also use an ansatz for 
the eq uatio n of state as suggested by IChevallier fc Polarskil 
l|200ll ) and lLinderl (|2003h (the CPL fit) 

w(z) =w + . (18) 

± ~~\~ z 

The results are shown in figure[9] We find, as expected, 
that it is impossible to catch the singularity in the equa- 
tion of state at z ~ 0.8 using an equation of state ansatz. 
Of course, one may try and improve upon this somewhat 
dismal picture by introducing fits with more free parame- 
ters. However, it is well known that the presence of more 
degrees of freedom in the fit leads to a larger degeneracy 
(between parameters) and hence t o larger errors of recon- 
struction (We ller fc Albrechtj r2002). In contrast to this ap- 
proach, when we reconstruct the equation of state using the 
smoothing scheme (which does not presuppose any particu- 
lar behavior of the equation of state), the Hubble parame- 
ter is reconstructed very accurately and hence the 'effective' 
equation of state for this model is also reconstructed well, as 
shown in figure From this figure we see a clear evidence 
of the singularity at z ~ 0.8. Thus to obtain maximum in- 
formation about the equation of state, especially in cases 
where the dark energy model is very different from the typi- 
cal quintessence-like models, it may be better to reconstruct 
the Hubble parameter or the dark energy density first. 

Therefore, we find that the smoothing scheme, which 
performs reasonably when reconstructing quintessence mod- 
els of dark energy models, can be also applied to models 



which show a departure from general relativistic behavior at 
late times. This section illustrates the fact that, in general, 
reconstructing H(z) and its derivatives such as the decelera- 
tion parameter q(z) may be less fraught with difficulty than 
a reconstruction of w(z), which, being an effective equation 
of state and not a fundamental physical quantity in some 
DE models, can often show peculiar properties. 



7 CONCLUSION 

This paper presents a new approach to analyzing supernova 
data and uses it to extract information about cosmological 
functions, such as the expansion rate of the universe H(z) 
and the equation of state of dark energy w(z). In this ap- 
proach, we deal with the data directly and do not rely on 
a parametric functional form for fitting any of the quanti- 
ties <Il(z), H (z) or w(z). Therefore, we expect the results 
obtained using this approach to be model independent. A 
Gaussian kernel is used to smooth the data and to calcu- 
late cosmological functions including H(z) and w(z). The 
smoothing scale used for the kernel is related to the number 
of supernovae, errors of observations and derived errors of 
the parameters by a simple formula, eq (|6]). For a given su- 
pernova distribution, the smoothing scale determines both 
the errors on the parameters and the bias of the results (see 
appendixfX]). A cannot be increased arbitrarily as this would 
diminish the reliability of the results. We use a value of A 
which gives results which have reasonably small bias as well 
as acceptable errors of H(z) for the SNAP quality data used 

3 Note, however, that most reconstruction methods including the 
present one may have problems in reproducing the rapidly oscil- 
lating e quation of state pred icted to arise in some models of dark 
energy jSahni fc Wanell2000h . 
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in our analysis (see section [3}- As can be seen from eq (|6]), 
when the data improves (i.e., the number of data points 
increases and/or measurement errors decrease), we expect 
that the same value of A would result in smaller errors on 
H{z). 

We demonstrate that this method is likely to work very 
well with future SNAP-like SNe data, especially in recon- 
structing the Hubble parameter, which encodes the expan- 
sion history of the universe. Moreover, our successful recon- 
struction of the Hubble parameter can also be used to dis- 
tinguish between cosmological models such as ACDM and 
evolving dark energy. The method can be further refined, 
if one wishes to reconstruct the cosmic equation of state to 
greater accuracy, by double smoothing the data- smoothing 
the Hubble parameter, after it has been derived from the 
smoothed luminosity distance, so as to reduce noise in w(z) 
(as in section [4}. The results obtained using the smooth- 
ing scheme compare favorably to results obtained by other 
methods of reconstruction. Another quantity which may be 
reconstructed to great accuracy is the look-back time of the 
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universe. 



An important result of this pa per is the di s cover y that 
the w-probe (originally proposed in lAlam et al.l l|2004 )) pro- 
vides us with an excellent diagnostic of dark energy. We 
summarize some of the attractive features of this diagnostic 
below. 

(a) The w-probe defined in Eqs. (|10[) and is ob- 
tained from the luminosity distance by means of a single dif- 
ferentiation. Therefore, it avoids the pitfalls of w(z) which 
is obtained from the luminosity distance through a double 
differentiation - see Eq.Q, a nd hence is usually accompa- 
nied b y large errors (see also iMaor. Brustein fc Steinhardtl 
(ill])). 

(b) The w-probe is robust to small uncertainties in 
the value of Qom- This attractive property allows us to get 
around observational uncertainties in the value of Qom cur- 
rently known to an accuracy of about 30%. Indeed, when 
marginalized over flom, the tu-probe can be used to great 
advantage to distinguish between ACDM and other dark en- 
ergy models for SNAP-quality data. 

We therefore conclude that the proposed reconstruction 
method by smoothing the supernova data appears to be suf- 
ficiently accurate and, when applied to SNAP-type observa- 
tions, should be able to distinguish between evolving dark 
energy models and a cosmological constant. 

The method proposed by us can also be used for other 
forms of data which deliver the luminosity (or angular size) 
distance. 
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APPENDIX A: SMOOTHING ERRORS AND 
BIAS 

In this section we explore the errors on the cosmological 
parameters due to the smoothing scheme, as also the bias 
which enters the results. 



Al Smoothing errors 

The smoothing scheme used in this paper is of the form : 

JV JV 

y(z) s =y G (z)+J2[y^i)-yG(zi)]S(z i ;A)/J2s(z i ;A) , (Al) 



The quantity being smoothed (in this case lnd^) is repre- 
sented by y, while y s represents the smoothed result. Let 
the errors in the data at any redshift Zi be given by a y (zi) 
and the errors in the guess model be a yG (zi). If we look at 
the second term on the right hand side of eq ()A1|I . we see 
that the errors on this term would be approximately given 
by the errors on y weighted down by the smoothing scale 
A and the number of data points N. Therefore the error on 
the smoothed result is: 



(z) 



,(*) + 



a 2 y {z) +o-y G {z) 
NA 



(A2) 



We now consider the errors for an iterative method. The 
first guess is an exact model, ACDM. Therefore the error on 
the result of the first iteration is simply 



<r vi (z) 



NA 



(A3) 



The next guess model is yi(z). Therefore the error on the 
result is 



o- 2 y2 (z) ~ [2 



(A4) 



J_\ o~l{z) 
NAJ NA ' 

From this we can show that the errors on the result for 
the M-th iteration is : 



M-l 

£ 



V 



iVA 



M + 



M(M-l) 
2NA 



°l(z) 
NA 



(A5) 



The second term on the right-hand side is small for a rea- 
sonable number of iterations, since N ~ 2000 and A > 0.01 
usually. Therefore we may approximate the errors on the log 
luminosity distance after M iterations for the guess model 



M 
Clndr 



w = vM ff ; dl ( z ) , 



(A6) 



where cr lndi (z) is the error for a simple smoothing scheme 
where the data is smoothed without using a guess model. 

A2 Smoothing Bias 

In any kind of a smoothing scheme for the luminosity dis- 
tance, some bias is introduced both in it and in derived 
quantities like H(z) and w(z). To illustrate the effect of this 
bias, we calculate it for the simplest Gaussian smoothing 
scheme for \ndi,(z) with the width A(z) <C 1 : 



lnd L (z) s = N(z) ^ \nd L (zi)exp 



2A 2 



N(z) 



M 

E 



exp 



2A 2 



(A7) 



(A8) 



where M is the total number of supernovae data points. The 
bias at each redshift (B(z) — lnd_t(z) s — lndi(z)) is the dif- 
ference between the smoothed lnd^^) and the exact value 
of lndi(z) : 



where the quantity S(zi; A) represents the smoothing func- 
tion with a scale A and ya(z) is the subtracted guess model. 



B(z) = N{z)^r(\nd L (zi)-\nd L (z)) exp 



2A 2 



.(A9) 
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z z 

Figure Al. Comparison of the reconstructed parameters obtained by using the smoothing method l)A7jl with the variable A(z) = 
0.2z/(l + z) 2 for the ACDM model (the dotted line in each panel) with the analytical biased result given by Eq. tATZl (dashed line). 
The solid line represents the fiducial ACDM model. Note the excellent agreement between the analytical and numerical results in the 
redshift range 0.1 < z < 1.1. 



Expanding lndi(aj) in terms of lnd,L(z) and its derivatives 
by Taylor expansion, we get: 

M r 



QndL(z))'(zi-z) + (\nd L (z))' 



, (Zi — z) 



x exp 



{Z-Zif 



2A 2 



,(A10) 



where the prime denotes the derivative with respect to z and 
we neglect higher derivatives. To see the effect of this bias at 
low and high redshifts where the number of supernovae on 
both sides of each z are not equal, we rewrite Eq. (|A10|) in 
another way. Let 5 be the spacing between two neighboring 
data points, so that z = mS. For m < M/2, we have: 

B(z) = 



N(z)J2 



(lndi(z)) 



„5 2 (i-m) 2 



exp 



M 

+N(z) 

i=2m+l 



(lnd L (z))'S(i - m) + {lnd L (z)) 



S 2 (i-m) 2 
2A 2 

„S 2 (i-m) 2 



x exp 



and for m > M/2: 
B(z) = 



N(z) 

i=2m-M+l 
2m-M 

+N(z) ^ 



(hid L (z)) 



„5 2 {i-mf 



exp 



S 2 (i-m) 2 
2A 2 



5 2 {i-rn) 2 



2A 2 



{\ad h {z))'5{i-m) + {\nd L {z))' 



,5 2 {i 



x exp 



8 2 {i-m) 2 
2A 2 



The first term in the above equations is the general bias 
of the method, while the second term is the bias arising 
due to an asymmetric number of data points around each 
supernova. For m — M/2, the number of data points is the 
same from both sides and we have: 



B(z)=N(z)J2 



{\nd L (z)) 



„6 2 (i-m) 2 



exp 



8 2 {i~mf 
2A 2 



In the continuous limit where x = i — m is assumed, we get: 

S 2 x 2 ' 



B(z) = N(z) I [(lnd L (z)) 
S 2 x 2 



s-2 2 
ll X 



J exp 



2A 2 



dx (All) 



N(z) = 



exp 



2A 2 



dx 



Therefore, the bias has the simple form 

K( (\nd L (z))"S 2 A 2 „ 
2^/A 2 =~( lnd£ W) ' 



(A12) 



This is a good analytical approximation for the bias at red- 
shifts in the middle range, where we do not encounter the 
problem of data asymmetry. To see the effect of this bias, 
let us assume that the real model is the standard ACDM, 
add the bias term to this model and then calculate the bi- 
ased H(z) and w(z). The result from this analytical calcu- 
lation can be compared to the result of smoothing the exact 
ACDM model using our method. The figure fATI simplv illus- 
trates that the results obtained using Gaussian smoothing 
and by the use of formula (|A12f) are in good agreement in 
the middle range of redshifts. However, we do not expect the 
formula (|A12|) to work properly at very low (z < 0.1) and 
high (z > 1) redshifts where the above mentioned asymme- 
try of points adds a further bias. 

Also, it appears that the smoothing bias has a tendency 
to decrease w(z) below its actual value in the middle range 
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of z. Thus, AC DM may appear to be a 'phantom' (w < — 1) 
if too large a smoothing scale is chosen. 



APPENDIX B: EXPLORING SMOOTHING 
WITH VARIABLE WIDTH A(Z) 

In order to deal with the problem of data asymmetry and 
paucity at low and high redshifts we may consider using 
a variable A(z). (i) Low z (z* C 1) : in this case, there 
are many more supernovae at z > z, than there are at 
2 < 2*. The error-bars are also small in the low redshift 
region. Therefore, a smaller value of A appears to be more 
appropriate at low z. (ii) High z (z» > 1) : in this case, 
there is considerably more data at z < z» than at z > 2*. 
However, at high z the errors are considerably larger than 
at low z, which suggests that in order to avoid a noisy result 
we must use a larger value of A in this region. In this sec- 
tion, we investigate two different functional forms of A(z) 
with the above properties and show how they result in the 
reconstruction of the equation of state. 



The results obtained using this method are shown in 
figure lB2l for SNAP data for the fiducial model w = —0.5. 
We find that this variable form of A(z) leads to a slight 
improvement of results at low redshifts by getting rid of the 
small bias which remains in the bootstrap iterative process. 
This improvement of the results is expected especially for the 
cosmological models whose equation of state at low redshift 
is very different as compared to the ACDM model, which is 
our initial guess model. 



Bl A(z) = A z/(1 + z) 2 

In section[2]we mentioned that, for \z— Zi\ <C 1, the exponent 
in Eq. © reduces to the form -(z - Zi) 2 /2A 2 (I + zf and 
the effective Gaussian smoothing scale becomes A(l + z). 
So if we use a variable A(z) — Aqz/(1 + z) 2 then the ef- 
fective Gaussian smoothing scale approaches a constant at 
large z and tends to a small value at small z. The results 
obtained using this method are shown in figure lBTl for SNAP 
data, using the model w = —0.5. We find that, the result 
for the Hubble parameter does not change much. However, 
the equation of state is somewhat better reconstructed, but 
noisier at low redshift because of the small width of smooth- 
ing. 



B2 tan-hyperbolic form of A(z) 

Tangent hyperbolic form for A(z) is another form of the 
variable A(z) which can simultaneously satisfy both the low 
and high z requirements. It has a small value at low redshifts 
and a bigger value at the higher redshifts. An additional im- 
portant property of this function is that it changes smoothly 
from low to high z, which translates into a smoother second 
derivative w(z) - see © - Q. 

A drawback of this method is that the tangent hyper- 
bolic function introduces a number of free parameters into 
the problem. However the role of these parameters can be 
understood as follows. The tangent hyperbolic function can 
be written in the general form 

A(z) = atanh^tf . (Bl) 

As we saw earlier, if A is held constant, then optimal results 
are obtained for Ao = 0.24 in © when we use bootstrap 
iterative process. We therefore determine a, b and c in (|B1|) 
so that A(z) w iA at z w 0, and A(z) « |A at z « 1.7; 
consequently 

A(z) = 0.36tanh^^y^ . (B2) 
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^ = 1-2^ 
Fiducial Model: w = -0.5 






Figure Bl. The smoothing scheme of equation ((2J is used with A(z) = X.2z/[1 + z) 2 to obtain smoothed H(z) and w(z) from 1000 
realizations of the SNAP dataset. The panel (a) represents the form of A(z) used, while panels (b) and (c) represent the reconstructed 
H(z) and w(z). The dashed line in panels (b) and (c) represents the fiducial 'metamorphosis' model with w = —0.5 while the solid lines 
represent the mean and la limits around it. The dotted line is ACDM. 



A = 0.36tanh^|±s 



Fiducial Model: w 
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Figure B2. The smoothing scheme of equation (f2j) is used with a tangent hyperbolic form of variable A(z) to obtain smoothed H(z) 
and ui(z) from 1000 realizations of the SNAP dataset. The panel (a) represents the form of A(z) used, while panels (b) and (c) represent 
the reconstructed H(z) and w(z). The dashed line in panels (b) and (c) represents the fiducial ui = —0.5 'metamorphosis' model while 
the solid lines represent the mean and la limits around it. The dotted line is ACDM. 
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